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Diffusion of electrons in a two-dimensional system with time-dependent random poten- 
tials is investigated numerically. In the absence of spin-orbit scattering, the conduc- 
tivity shows universal weak localization correction. In the presence of it, however, the 
correction to the conductivity weakly depends on the strength of disorder, and becomes 
vanishingly small close to the metal-insulator transition point. 

1. Introduction 

Transport properties of the two-dimensional disordered electron systems have attracted much at- 
tention in recent years. To understand the transport properties of a random system, the concept 
of quantum interference plays an important role. Effects of quantum interference are suppressed 
in the presence of the inelastic scattering, for instance, due to the electron-phonon interaction, 
due to the Coulomb interaction of conduction electrons, or due to the motion of a single impurity 
atom. The purpose of this paper is to calculate numerically the conductivity of the two-dimensional 
systems with time-dependent random potentials with and without spin-orbit scatterings, and to 
demonstrate dephasing by the time-dependent potential. 

The two-dimensional system with spin-orbit interactions belongs to the symplectic universality 
class, and exhibits the Anderson transition. On the other hand, it belongs to the orthogonal 
universality class without spin-orbit scatterings. The study of the effect of dynamical potential 
on the conductivity in the critical regime as well as the weakly localized regime gives how the 
dephasing modifies the quantum interference effects in the above two universality classes. 

2. Model and Method 

We consider the tight-binding Hamiltonian with time-dependent potential £{{t) on the two-dimensional 
square lattice: 



where Cj a (Ci j(7 ) denotes a creation (annihilation) operator of an electron at the site i with spin a 
and are the nearest neighboring sites. Here x(y) denotes the unit vector in the x(y)— direction. 
All the length-scales are measured in units of the lattice constant ag. The strength of the spin-orbit 





with the transfer integral in the Ando model Q 





interaction is characterized by the parameter S = V2/V with V 




In order to take into account the effect of moving potential, we assume that the site-potentials 
take the form: 

e f (i) = e r (O)V^cos[a;i + 0(r)], (3) 

where uj is the frequency. Effects of scattering from impurities are introduced through randomness 
of site energy ey(0) and phase 6{f) at t = distributed uniformly in the regions, |£f(0)| < W/2 and 
\0(r)\ < 7r. We consider the adiabatic case u <C V/K, where impurities move slower than electrons. 
The method of the calculations is essentially the same as those in Refs. |2Lj| where we adopt the 
method based on the fourth-order decomposition of exponential operators Qj. 

The quantity we observe is the second moment of the wave packet {r 2 (t)) c defined by 

(r 2 (t)) c =(\f(t)\ 2 )-(r(t)) 2 (4) 

with 

(r 2 (t)} = /dOr d-1 drr 2 |^(f,t)| 2 , (5) 

where ijj(r, t) denotes the wave function at time t, and d the dimensionality of the system. If the 
wave function extends throughout the whole system, the second moment is expected to grow in 
proportion to time t as (r 2 ) c = 2dDt. Here, the coefficient D denotes the diffusion coefficient. In 
contrast, if the wave function is localized, it is clear that the second moment remains finite in the 
limit t —> oo|J. In the metallic region the Einstein relation a = e 2 p(Ep)D relates the conductivity 
a to the diffusion constant, where p(Ep) is the density of state at Fermi energy Ep. 

3. Numerical Results 

We have calculated the second moment of the wave packet (r 2 ) c at various random potential 
strength. The size of the systems are 500 by 500 for the energies Ep/V = —1. We have carried 
out an exact diagonalization for the 20 by 20 subsystem at the center of the system and taken the 
eigenfunction of the subsystem whose eigenvalue is closest to the given energy Ep as the initial 
wave packet. By this procedure we can simulate the diffusion of the wave packet whose energy is 
approximately equal to Ep. 

Fig. 1 gives an example of the calculated second moment for uj = W~ 2 V/h. The second moment 
is proportional to the time in a wide range up to 2xW 3 h/V. We can evaluate the diffusion constant 
of this system with least square fit to these data. In the actual simulation, the quantities (r 2 ) c 
are averaged over at least four samples of random potential distribution. The density of states 
p(Ep) in the Einstein relation is evaluated from the direct diagonalization or numerical Green's 
function method for uj = 0. We have obtained that p{Ep)alV = 0.281 and 0.251 for W/V = 2^2 
and 4, respectively in the orthogonal system 5 = 0, and p{Ep)a 2 ] V = 0.304,0.267 and 0.222 for 
W/V = 2\/2, 4 and 4v2, respectively in an example of the symplectic system 5 = 0.5. The static 
conductivity obtained by the present method coincides with that by Landauer formula 

In Fig. 2 we show examples of the conductivity in the presence (open circles) and absence 
(filled circles) of spin-orbit interaction as the function of In a) where uj = lo j {V jfi). In the absence 
of spin-orbit interaction the conductivity linearly depends on In a; as 

a = Alnu + <r(0). (6) 

The coefficients ^4(~ 0.22) calculated by the least square fit to the averaged conductivity, which are 
shown by the dotted lines, are almost the same among these different disorder cases, although the 
conductivity cr(0) for tu = is significantly different from each other. For 5 = 0.5 the conductivity 
decreases linearly, takes a minimum on lnu; ~ —2, and then increases with In a). For W/V = 2^/2 



and 4, the absolute value of the slopes on lnu) < —2 is almost half of orthogonal case. The slope 
on lncD < —2 represented by dashed lines slightly decreases with the increase of the disorder W/V, 
and then it becomes almost flat for W/V = 4^2. Note that, W/V = 4^2 = 5.66 . . . is near the 
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Fig. 1. An example of the second moment as the 
function of time for to = 10~ 2 V/h and E-p/V = 
— 1. The dashed line is that in the absence of 
spin-orbit interactions (5 = 0) and the solid lines 
in their presence (S = 0.5). The evaluated diffu- 
sion constants with least square fit to these data 
are Dh/a 2 V = 3.49 and 1.67 for W/V = 2^2 
and 4 in S = 0, and Dh/a 2 V = 4.06,2.11 and 
1.06 for W/V = 2^/2,4 and 4^2 in S = 0.5, 
respectively. 
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Fig. 2. The calculated conductivity of the two- 
dimensional system with time dependent im- 
purity potential for W/V = 2^/2,4 and 4\/2. 
The dotted lines indicate in orthogonal system 
the corrections 5a = 0.217^- In u), and 5a = 

2 

0.214^- In u, and the dashed lines in symplec- 
tic system the corrections on lnu) < — 2 5a = 
-0.138^ In u, 5a = -0.100^ In cu, and 5a = 

2 

—0.017^- In cj, respectively. 



critical point of metal-insulator transition W c /V 
for E-p = OQ. 



5.74 estimated by the finite-size scaling method 



4. Discussion and Conclusion 

To interpret the above universal correction due to frequency uj, we recall the weak localization 
correction, 



5a 



■^-ln^ 

27T 2 h T 



(7) 



where is the dephasing time, and r the elastic scattering time. The u dependence of T& in 
the kicked rotator problem || as well as that in the one-dimensional disordered system || in the 
presence of noise is estimated to be 

t> ~ uj- 2 / 3 . (8) 

Since the qualitative argument leading to this dependence is independent of the space dimension, 
the 2D weak localization correction 5a ||] is estimated to be 



5a 



2 e 2 uj 
3ir h ujq 



(9) 



The pre-factor 2/37r = 0.212... of the In a; term agrees with the numerical calculation shown in Fig. 
2, which is universal and independent of the potential strength W. Note that this (^-dependence is 
also observed in the case of low energy phonon scattering || . 

The results with the spin-orbit interaction is consistent with the weak localization correction 
too. For T0 <C tso, where tso is the dephasing time due to spin-orbit interaction, the conductivity 
increases with to as in the orthogonal case, while the correction to the conductivity is half with 
opposite sign for 3> tso- 

In conclusion, we have analyzed the weak localization effect in the two-dimensional system with 
time-dependent random potentials using the equation-of-motion method. It has been shown nu- 
merically that the weak localization correction to the conductivity due to the fluctuating potentials 
takes universal value independent of the random potential strength W in orthogonal system. The 
symplectic system is also investigated and it is demonstrated that the conductivity becomes almost 
independent of the frequency near the critical point of metal-insulator transition. Such behavior is 
observed in the recent experiments with high mobility Si MOS |R]], though the mechanism leading 
to the metal-insulator transition is different. Our results will open a new way to incorporate the 
dephasing mechanism in numerical simulations. 

Acknowledgments 

Numerical calculations were performed on FACOM VPP500 in Supercomputer Center at RIKEN 
and ISSP. 

References 

[1] T. Ando, Phys. Rev. B40 (1989) 5325. 

[2] T. Nakanishi, T. Ohtsuki, and T. Kawarabayashi, J. Phys Soc. Jpn. 66 (1997) 949. 

[3] T. Kawarabayashi and T. Ohtsuki, Phys. Rev. B51 (1995) 10897; B53 (1996) 6975. 

[4] M. Suzuki, Phys. Lett. A146 (1990) 319; A165 (1992) 387; J. Math. Phys. 32 (1991) 400; J. 
Phys Soc. Jpn. 61 (1992) 3015; Proc. Jpn Acad. 69, Ser. B (1993) 161; Commun. Math. Phys. 
163 (1994) 491. 

[5] S. Fishman and D. L. Shepelyansky, Europhys. Lett. 16 (1991) 643. 
[6] F. Borgonovi, and D. L. Shepelyansky, Phys. Rev. E51 (1995) 1026. 
[7] G. Bergmann, Phys. Rep. 107 (1984) 1. 

[8] S. Hikami, A. Larkin and Y. Nagaoka: Prog. Theor. Phys. 63 (1980) 707. 

[9] In the case of low energy phonon scattering, the dephasing time is estimated to be = 
(3T/2mE F ) 1/3 v- 2 / 3 , where v is the phonon velocity. Assuming uocwwe obtain ~ u) ' . 
See, M. J. Stephen, Phys. Rev. B36 (1987) 5663. A. A. Golubentsev, Zh. Eksp. Teor. Fiz. 86 
(1984) 47. [Sov. Phys. JETP 59 (1984) 26.] 

[10] S. V. Kravchenko, W. E. Mason, G. E. Bowker, J. E. Furneaux, V. M. Pudalov, and M. D'lorio, 
Phys. Rev. B51 (1995) 7038. 



